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Abstract 

The search for a general means to control the coupling between electricity and magnetism 
has intrigued scientists since 0rsted's discovery of electromagnetism in the early 19 th century. 
While tremendous success has been achieved to date in creating both single phase and composite 
magnetoelectric materials, the quintessential electric-field control of magnetism remains elusive. 
In this work, we demonstrate a linear magnetoelectric effect which arises from a novel carrier- 
mediated mechanism, and is a universal feature of the interface between a dielectric and a spin- 
polarized metal. Using first-principles density functional calculations, we illustrate this effect at 
the SrRuOs/SrTiOs interface and describe its origin. To formally quantify the magnetic response 
of such an interface to an applied electric field, we introduce and define the concept of spin ca- 
pacitance. In addition to its magnetoelectric and spin capacitive behavior, the interface displays 
a spatial coexistence of magnetism and dielectric polarization suggesting a route to a new type of 
inter facial multiferroic. 
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The linear magnetoelectric effect is the first-order magnetic response of a system to an 
applied electric field, or equivalently the electrical polarization induced by an applied 
magnetic field 1,2 : 
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where a. is the magnetoelectric tensor, which is non-zero in the absence of time-reversal 
and space-inversion symmetries. Although the effect has long been recognized both 
theoretically 3 and experimentally 4 it has traditionally been of only academic interest be- 
cause of its inherent weakness in intrinsic single phase materials. The last few years, 
however, have seen a tremendous revival of activity in magnetoelectrics 2 , motivated in 
large part by the entirely new device paradigms that would be enabled by electric-field- 
control of magnetism 5,6 . Progress has been fueled by the growth of novel single-phase^ 
and composite 9,10 multiferroic materials, which show simultaneous magnetic and ferro- 
electric ordering and have the potential for significantly enhanced magnetoelectric re- 



currently two mechanisms for magnetoelectric coupling are well established. In tra- 
ditional single-phase magnetoelectrics, of which Cr2C>3 is the prototype, an intrinsic mag- 
netoelectric coupling occurs 3 . An electric field both shifts the positions of the magnetic 
cations relative to the anions - changing both the dipolar and exchange contributions 
to the magnetic interactions - and modifies the electronic wavefunctions, which further 
changes the magnetic coupling. This behavior, and its inverse when a magnetic field 
is applied, are enhanced in materials with strong spin-orbit coupling (and hence strong 
coupling of the magnetic order to the lattice) and with large piezo- or electro-strictive re- 
sponse. In composites of magnetostrictive or piezomagnetic materials combined with 
electrostrictive or piezoelectrics, an extrinsic magnetoelectric coupling is mediated by 
strain at the interface 9 ^. Here an electric field causes strain in the electrical compo- 
nent which is mechanically transferred to the magnetic component where it changes the 
magnetization (and vice versa) 44 . This class of structures appears very promising for 
technological applicability, in part because the underlying working principles are essen- 
tially based on the bulk properties of the constituents. Very accurate modeling has thus 
been achieved at the level of semiclassical theories 9 , which has allowed for efficient opti- 



sponses. 



2 



m an e 
Ref. 131 



mization of geometry and material properties. 

Recently, it has been pointed out that thin film growth of heterostructures can be ex- 
ploited to create structures which break inversion symmetry even when the parent bulk 
materials are centrosymmetric 12 . In fact at the interface between any two unlike mate- 
rials the space inversion symmetry is intrinsically broken. Therefore, provided that one 
of the two components breaks time-reversal symmetry (which occurs in the case of a 
ferromagnet), we should expect a linear magnetoelectric coupling to occur, even in the 
absence of a piezo- or electro-strictive material (Table fl). An intriguing recent observa- 
tion of electric-field tunable magnetization at the surface of a ferromagnetic electrode 
ectro-chemical cell 13 may indeed be a manifestation of such an interfacial effect, 
suggested an important role played by accumulation/ depletion of carriers at the 
ferromagnetic surface, but a detailed microscopic mechanism could not be determined. 
A simpler and more convenient architecture for developing a quantitative and general 
model is a capacitor geometry, in which one of the two components of the bilayer is a 
metal and the other is an insulating dielectric, and a second electrode is deposited on top 
of the insulating film to provide a source of external bias; we study this geometry here. 

The magnetoelectric response of such an interface is intimately related to the local de- 
tails of chemistry, bonding, structure and electronics, which often strongly depart from 
those of single-phase compounds; to account for the subtle balance of these often com- 
peting factors requires a fully quantum-mechanical treatment. Density functional theory 
has proven to be a very effective tool to deal with the challenging problem of complex 
oxide interfaces; however, the lack of a rigorous methodological framework to treat finite 
electric fields in metal-insulator heterostructures has until recently thwarted attempts to 
access dielectric and magnetoelectric responses of such systems. Recent advances in fi- 
nite electric field methods^^, which we extended in this work to include the magnetic 
degree of freedom, enables such calculations for the first time. 

Here we demonstrate computationally a novel linear magnetoelectric effect which is 
carrier-mediated in origin. The effect occurs at the interface between a non-magnetic, non- 
polar dielectric and a metal with spin-polarized carriers at the Fermi level. Using first- 
principles density functional theory, we find that 

(i) an electric field induces a linear change in magnetization in the metal at the inter- 
face, with the sign of the change determined by the orientation of the applied field, 



(ii) the magnetic response is mediated by the accumulation of spin-polarized charge 
within the few atomic layers closest to the interface; this charge is stored capaci- 
tively leading to a spintronic analogue to a traditional charge capacitor, and 

(iii) as a result of the magnetoelectric response, magnetism and dielectric polarization 
coexist in the interfacial region, suggesting a route to a new type of interfacial multi- 
ferroic. 

We choose SrRu03/SrTi03 heterostructures as our representative system, motivated 
by the favorable dielectric properties of SrTi03 (STO) 16 and the widespread use of thin 
film ferromagnetic SrRu03 (SRO) as oxide electrodes 17 ; furthermore, the effects of interfa- 
cial carrier modulation in thin films of both materials have been intensively investigated 
as a possible route to the realization of novel field-effect devices 1 ^*^. We note that our 
conclusions are not specific to this system but should apply to any dielectric /magnetic 
metal interface. Since such interfaces are similar to those already in place for magnetic 
tunnel junctions 20 , rapid experimental verification of this effect should be well within 
reach. 

Our calculations are performed within the local spin-density approximation (LSD A) 
of density functional theory and the projector-augmented wave (PAW) method 45 , with 
a planewave basis cutoff energy of 40 Ry. Our model heterostructure consists of seven 
unit cells of SrRu03 periodically alternating with seven unit cells of SrTi03 (Figure [1]) 
along the (100) direction. The in-plane lattice parameter was set to the theoretical SrTi03 
equilibrium lattice constant (3.85 A) to simulate epitaxial growth on a SrTi03 substrate. 
The two-dimensional Brillouin zone of the resulting tetragonal cell was sampled with six 
special fc-points, corresponding to a 6 x 6 Monkhorst-Pack mesh 22 . The ion positions and 
out-of-plane lattice constants were fully relaxed until the forces on the ions and the out- 
of-plane stress were less than 1 meV A -1 using an extension of the Car and Parrinello 
method for metallic systems 2 ^* 2 ^ 2 ^. To study dielectric and magnetoelectric responses 
we used a spin-polarized extension of our recently developed method for applying fi- 
nite electric fields to metal-insulator heterostructures^^ 5 -. With these parameters we ob- 
tain structural, dielectric and electronic properties for bulk SrTi03 and SrRu03 that are 
consistent with earlier theoretical studie s 26 f27f28f29 . In particular, our calculated dielectric 
constant for cubic SrTi03 at this lattice constant is e = 490, in excellent agreement with 
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Ref . 26, and our calculated magnetic moment for ferromagnetic SrRuC>3 epitaxially con- 



strained to the LDA equilibrium lattice constant of SrTi03 (less than 0.2% strain) is 1.03 
He per Ru atom, consistent with Refs. 



27 



28 



29. 



To investigate the response of the heterostructure to an applied electric field, we ap- 
ply an external bias of AV =27.8 mV across the capacitor plates, and allow the ions to 
relax to their equilibrium positions. The resulting planar and macroscopically averaged 
change in magnetization is shown in Figure [0 The overall induced magnetic moment 
is localized at the interfaces, and amounts to 2.5 x 10~ 3 jib per surface unit cell, corre- 
sponding to a surface spin density of 0.27 ^C/cm 2 . The accumulation of spin is exactly 
equal in magnitude and opposite in sign (within the numerical accuracy of our method) 
at the left and right electrode so that the overall induced magnetic moment of the het- 
erostructure is zero, consistent with the symmetry of the system. In the same Figured] we 
show the planar average of the induced spin density without macroscopic averaging. It 
is clear that the dominant contribution to the induced magnetization is accumulation of 
spin on the interfacial RuC>2 layer; interestingly, this effect is partially compensated by a 
smaller opposite induced spin density in the adjacent RuC>2 layer, somewhat reminiscent 
of Kondo behavior 30 . Conversely, the induced spin on the SrTiC>3 side of the interface, 
which is provided by the exponentially vanishing tails of the metal-induced gap states, 
is small. 

Since the magnetoelectric tensor, a, is an intrinsic bulk quantity which relates bulk 
induced magnetization to applied electric field, it is not the appropriate parameter for 
describing the interfacial effects observed here. Instead, to characterize the magneto- 
electric performance of the device, we start by introducing the notion of spin-capacitance 
density per unit area, C s which we define as 

C s = (3) 

where o~ s is the amount of spin polarization per unit area induced by the voltage V. This 
definition parallels that of the charge capacitance density per unit area, C = y, where a 
is the surface density of free charge stored at the electrode. Capacitance densities, how- 
ever, depend both on the material properties and on the capacitor geometry, in particu- 
lar, on the thickness of the dielectric film. Therefore, in order to obtain a parameter that 
characterizes the fundamental properties of the interface response, we further define the 



dimensionless parameter r\, which is the ratio of spin capacitance to charge capacticance, 
i.e. r] = C s /C. By definition, rj is also equal to the ratio of the induced spin with respect to 
the free charge accumulated at the capacitor plates, and is clearly zero for non-magnetic 
electrodes; this interpretation provides a very intuitive picture of the interface response 
to a field. The spin capacitance of our structure is C s =96.0 £F/ um 2 , to be compared with 
a charge capacitance of C°=258 fF/ jim 2 , yielding a value of //°=0.37. 

In order to investigate the origin of our calculated linear magnetoelectric effect, we 
first compare this value of C° (calculated here by using spin-polarized SRO electrodes) 
with the value obtained in Ref. La for the same system calculated within the non-spin- 
polarized local density approximation. Strikingly, the two capacitance values are iden- 
tical within numerical accuracy, indicating that the presence (or absence) of magnetic 
ordering in the electrode does not influence the dielectric response of the capacitor. The 
static dielectric response of an insulator can be decomposed into an electronic and an 
ionic contribution, where the latter is related to the frequency coi and the dipolar activity 
of the zone-center phonons: 

Here Coo is the average frozen-ion permittivity of the supercell, O is the volume and Z, 
are the "mode effective charges" obtained by summing the ionic Born effective charges 
weighted by their relative contributions to the normal mode eigenstates. We found no ap- 
preciable differences in either €<*> or the Born effective charges between the non-magnetic 
and the spin-polarized calculations; the fact that the final value of C° is also unaffected 
strongly suggests that the frequencies of the dipolar-active modes are themselves insen- 
sitive to magnetism. Therefore, coupling between magnetism and phonons, which has 
been recognized as a viable route for magnetoelectric switching in perovskite materials 31 , 
is not the driving force for the magnetoelectric response observed here, and a different 
effect must be responsible. 

To elucidate this behavior, we repeat our electric field calculations with the ions frozen 
in their initial centrosymmetric positions (Figure |2]); this isolates the electronic response 
of the system and corresponds experimentally to the high-frequency limit. Since by keep- 
ing the ions fixed we completely remove any structural effect on magnetism, any induced 
magnetization must be of intrinsic electronic origin, and can only be ascribed to the ca- 
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pacitive accumulation of spin-polarized carriers at the interface. In this case we obtain an 
induced magnetic moment at each interface of 1.80 x 10 -4 y.^ P er un it cell, corresponding 
to a spin-polarized charge density of 0.019 ^C/cm 2 , and a high frequency spin capaci- 
tance of Cf= 6.99 fF/^m 2 . All of these values are an order of magnitude lower than the 
corresponding quantities in the static case. Indeed the substantially lower induced mag- 
netization and spin capacitance parallel an order of magnitude reduction in the charge 
capacitance, C°°=20.3 fF/ jim 2 , which in turn derives from the suppressed dielectric re- 
sponse of the system when the ions are not allowed to relax. As a consequence, the 
corresponding value of t/°°=0.34 is very close to that obtained in the static case. The simi- 
larity in the values of rj calculated in the static and high-frequency limits indicates that the 
mechanisms leading to the screening of polar phonons in the static regime are the same 
as those screening electronic bound charges in the high-frequency regime, where changes 
in chemical bonding and /or structure are not possible. Therefore we conclude that our 
system provides the first example of a carrier-mediated magnetoelectric effect, which results 
entirely from the capacitive accumulation of spin-polarized carriers at the interfaces. 

The origin of the carrier-mediated magnetoelectric effect is shown schematically in 
Figure |3j On application of an external field, free carriers accumulate at the capacitor 
plates, which are partially screened by the dielectric polarization of the STO film. In 
the half -metallic limit all displaced electrons are spin-polarized in the same direction (up 
in the figure); in the present case there is a partial cancellation between spin-up and 
spin-down carriers that reflects the incomplete spin polarization at the Fermi level of 
the interfacial SRO layer. This process accumulates up-spin magnetization adjacent to 
the positively charged electrode, leaving behind an absence of up-spin magnetization (or 
equivalently down-spin polarized holes) at the negative plate. 

One possible mechanism for detecting the induced magnetization is the magneto- 
optical Kerr effect. Indeed, we note that the recent report of enhanced Kerr response 
in a 2 nm thick FePt film when an electric field is applied across an electrochemical cell 
with FePt electrodes 13 is also likely a manifestation of the carrier-mediated magnetoelec- 
tric effect that we describe here. In Ref. [13| a Kerr rotation of ~1° was obtained for an 
applied potential of 600 mV, which generated 0.015 electrons per unit cell; in our case 
an equivalent potential would induce a spin polarization of 0.05 spin-polarized electrons 
per unit cell, which should therefore be readily detectable. 



7 



Substitution of our field-polarized dielectric by a ferroelectric with a spontaneous elec- 
trical polarization would yield a correspondingly larger spin density and interfacial mag- 
netoelectric effect as well as an added degree of freedom from the switchable polariza- 
tion. The parameter n has a particularly transparent meaning in this case, relating the 
accumulated spin at the interface to the bulk polarization of the ferroelectric, P: Ne- 
glecting the effect of depolarizing fields (which is a good approximation in the limit of 
a sufficiently thick ferroelectric film), the free charge accumulated at the electrode ex- 
actly compensates the surface polarization charge induced by the spontaneous electric 
displacement. Therefore, r\ = cr s /P, and the induced spin at the interface can be read- 
ily predicted by the knowledge of n and P, under the assumption that the behavior of 
the free carriers is linear in P. For example, from the numbers calculated here we pre- 
dict a s =12 jWC/cm 2 for a SrRu03/BaTi03 interface (Pgro=32 uC/cm 2 ); this corresponds 
to a fully switchable magnetic moment of 0.11 ^g. To substantiate these arguments we 
repeated the calculation by replacing STO with BTO 46 . We restricted our analysis to the 
high-frequency case to avoid complications related to the proper treatment of BTO lat- 
tice instabilities; by analogy with the SRO/STO case, it is unlikely that the fully relaxed 
static behavior will significantly differ. The resulting 7/^ TO = 0.37 is indeed very close 
to the SRO/STO value, confirming the validity of our arguments. In addition, such fer- 
roelectric/ferromagnetic metal heterostructures would exhibit a new type of interfacial 
multiferroism, resulting from the simultaneous small magnetization of the ferroelectric by 
the magnetic metal reported in this work, and the polarization of the metallic magnet by 
the depolarizing field of the ferroelectric described previously in Ref . [3. We note that 
the interfacial multiferroism results from the attempt by the metallic electrode to screen 
the capacitive charge at the interface, and is mediated by the capacitive accumlation of 
spin-polarized carriers at the interfacial layers. The behavior is therefore distinct from 
that in composite multiferroics, which combine ferro- or ferri-magnetic insulators with 



ferroelectrics (see for example Refs. 9. lO Jill) . 

This class of heterostructures was recently investigated by Duan et al.—, who calcu- 
lated changes in magnetism induced by ferroelectric displacements at the interface be- 
tween ferroelectric BaTi03 and ferromagnetic Fe. According to Ref. 32, the positively and 
negatively polarized faces of the BaTi03 film have magnetic moments at the interface 
with the Fe electrodes which differ by almost 0.3 }Ib per surface unit cell; the differences 
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are ascribed to changes in interracial chemical bonding and hybridization. While the 



magnitude of this va 
in moment in Ref . 



ue is consistent with our estimates, surprisingly most of the change 
321 resulted from changes in spin polarization of the Ti d states, which 
we find to be minimally affected by the electric field. We note, however, that such changes 
in hybridization are strongly dependent on the position of Ti d states with respect to the 
Fermi level and can be affected by the well-known underestimation of the fundamen- 
tal insulating gap that plagues most commonly used DFT functionals. Conversely, the 
carrier-mediated accumulation of spin-polarized carriers investigated here is a general 
characteristic of all ferromagnetic electrodes in contact with insulators, regardless of the 
details of the bonding; since its magnitude is similar to hybridization effects, it should 
not be neglected. 

Finally, we mention some possible related applications that could exploit the predicted 
carrier-mediated magnetoelectric behavior. With the ability to capacitively store spin- 
polarized charge, we see that this heterostructure provides the spintronic analogue of a 
thin film charge capacitor, and could therefore find application in spin-polarized exten- 
sions to regular capacitor applications, such as a filter of spin-polarized direct current, or 
in spin logic circuits. In addition, since the spin capacitance mechanism only allows cur- 
rent to flow when the leads are magnetized in the same direction, the capacitor could be 
used as a sensitive magnetoresistive detector with the desirable feature that the electrons 
never cross interfaces and so are less susceptible to scattering of their spin orientation. 
Since the change in total magnetization is small this material combination is not immedi- 
ately transferable to electric field tunable inductors. The magnitude of the effect would 
be enhanced, however, by the use of materials with higher spin polarization at the Fermi 
surface, and lower overall magnetization; in particular half-metallic antiferromagnetic 
electrodes 33 with zero net magnetization would give an infinite relative change in net 
magnetization with electric field. We note also that the calculated LSDA dielectric con- 
stant of SrTiC>3 is 490, whereas experimental values can be tens of thousands, and the 
induced magnetization will scale accordingly. 

To summarize, we have used first-principles density functional calculations combined 
with finite electric field methods to demonstrate a linear magnetoelectric response in a 
complex oxide heterostructure. This new carrier-mediated mechanism adds an addi- 
tional degree of freedom in the design and functionalization of magnetoelectric multi- 



9 



ferroic materials, and may guide the way forward in designing electric-field switchable 
magnetic devices. With growing interest in superlattice multiferroic composites, the phe- 
nomenon we describe may elucidate experimentally observed effects that to date are at- 
tributed to unknown interfacial responses to applied electric (and magnetic) fields 34 . Our 
results also expand the growing body of literature demonstrating the novel functionality 
that can be achieved at oxide interfaces 1 ^ 13 ^ 3 ^ 13 ^ 13 ^ 3 ^ 4 ^ 41 and suggest an additional route 
to novel oxide-based interfacial devices 42 . 
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Although not strictly a linear magnetoelectric effect, a related behavior in single-phase multi- 
ferroics, in which electric field re-orientation of ferroelectric polarization reorients the magnetic 
easy axis or plane through the elastic coupling 43 and vice versa^, has recently received consid- 
erable attention. 

The PAW pseudopotentials were generated in the 4s 2 4p 6 5s 2 configuration for Sr, 3s 2 3p 6 4s 2 3d 2 
for Ti, 2s 2 2p 4 for O and 4s 2 4p 6 4d 7 5s 1 for Ru. 

Six Sr atoms within the insulating layer were replaced with Ba, which corresponds to assuming 
a T1O2 (SrO) termination of the BTO (SRO) lattice. The in-plane lattice parameter was kept 
fixed to the calculated bulk STO value, while both out-of-plane lattice parameter and atomic 
positions were fully relaxed while imposing mirror symmetry with respect to a plane parallel 
to the interface. 
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FIG. 1: (color online) Calculated magnetization induced by an external voltage of 27.8 mV in a 
nanocapacitor consisting of 7 layers of SrTiOa alternating with 7 layers of ferromagnetic, metal- 
lic SrRuC>3. The light (red) and bold (black) lines show respectively the planar averaged, and 
macroscopically and planar averaged induced magnetization profiles. 
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FIG. 2: (color online) Calculated high-frequency planar averaged magnetization induced by an 
external voltage of 27.8 mV. The light (red) and bold (black) lines show respectively the planar 
averaged, and macroscopically and planar averaged profiles. Note that the y-axis scale is exactly 
one order of magnitude smaller than in the static case (Fig. 1). The strong change in local spin 
density on the interfacial Ru atoms is evident, along with a small local magnetic moment on the 
interfacial Ti atoms. 
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FIG. 3: (color online) Schematic of the mechanism causing the induced magnetization at a fer- 
romagnetic (FM) metal - non-polar (NP) dielectric interface. The inset density of states cartoons 
show the changes to the occupation of half-metallic densities of states in the electrodes when the 
current flows in response to the applied electric field. The accumulation of up-spin electrons ad- 
jacent to the positive face of the dielectric, and their depletion from the negative face, leads to the 
net magnetization shown in Figs. [T] and |2l 
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m a cm ofi (-• 
llLdiil Lt; L1L. 


r\ r\r\ — tyi a cm ta-H n 


lILctiil Lt; Lit. 


llLcLiil Lt; L1C 




non-polar 


polar 


polar 


interface 


time-reversal symmetric? 


no 


yes 


no 


no 


space-inversion symmetric? 


yes 


no 


no 


no 


linear ME effect allowed? 


no 


no 


yes 


yes 



TABLE I: The linear magneto-electric (ME) effect is only allowed in the absence of both time- 
reversal and space-inversion symmetries. For a single phase material the co-existence of electrical 
polarization and magnetic ordering are required for linear magneto-electricity to occur; in the 
presence of an interface which breaks space-inversion symmetry the existence of magnetic order- 
ing is sufficient. 
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